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Abstract. The integral equations for the correlation functions of an inhomogeneous 
fluid mixture are derived using a functional Taylor expansion of the free energy around 
an inhomogeneous equilibrium distribution. The system of equations is closed by the 
introduction of a reference functional for the correlations beyond second order in the 
density difference from the equilibrium distribution. Explicit expressions are obtained 
for energies required to insert particles of the fluid mixture into the inhomogeneous 
system. The approach is illustrated by the determination of the equation of state of 
a simple, truncated Lennard- Jones fluid and the analysis of the behavior of this fluid 
near a hard wall. The wall-fluid integral equation exhibits complete drying and the 
corresponding coexisting densities are in good agreement with those obtained from 
the standard (Maxwell) construction applied to the bulk fluid. Self-consistency of the 
approach is examined by analyzing the virial/compressibility routes to the equation 
of state and the Gibbs-Duhem relation for the bulk fluid, and the contact density 
sum rule and the Gibbs adsorption equation for the hard wall problem. For the 
bulk fluid, we find good self-consistency for stable states outside the critical region. 
For the hard wall problem, the Gibbs adsorption equation is fulfilled very well near 
phase coexistence where the adsorption is large. For the contact density sum rule, 
we find deviations of up to 20% in the ratio of the contact densities predicted by the 
present method and predicted by the sum rule. These deviations are largely due to 
a slight disagreement between the coexisting density for the gas phase obtained from 
the Maxwell construction and from complete drying at the hard wall. 
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1. Introduction 

Integral equations have become a popular and increasingly accurate tool in describing 
the thermodynamics and the bulk structure of one-component simple liquids and, to a 
smaller degree of accuracy, of fluid mixtures. In the early, formal work (see references 
PU El Ej and references therein) the structural ingredients (pair correlation function, 
direct correlation function and bridge function) have been defined in terms of graphical 
expansions and two exact equations relating these three functions have been established. 
The unknown third equation which is necessary for a closed system of equations is 
usually referred to as the closure relation. A number of approximations have been 
developed in the past for the closure relation. Some of the more succesful ones for the 
structure and thermodynamics of a one-component liquid are the optimized random 
phase approximation (ORPA) the hierarchical reference theory (HRT) jHj, the 
self-consistent Ornstein-Zernike approximation (SCOZA) j^J, the Martynov-Sarkisov 
closure p[j and the reference hypernetted chain equations (RHNC) 8\. For a general 
introduction into the subject treating some of these and other common closures, see 
reference jO]. The extension of integral equation theories to a general inhomogeneous 
system defined by an external potential (a wall, say) is possible by incorporating the 
source of inhomogeneity as a second component in the dilute limit (the wall particle) 
and treating the originally inhomogeneous problem as a bulk (homogeneous) problem 
in the two-component mixture. This generalization fails for all of the above closures 
whenever the phenomenon of wetting in the presence of the external (wall) potential is 
involved. 

On the one hand this failure is linked to the difficulties of defining the free energy 
and chemical potential accurately and consistently in the two-component (wall-fluid) 
mixture. Or, whenever a free energy approximation underlying a certain closure relation 
can be formulated, its form does not allow for wetting. For example, it can be proved 
that in the hard wall case the HNC and RHNC closures miss the phenomenon of complete 
drying for the above reason (for HNC, the proof is given in reference ^Hj and it can be 
easily extended to RHNC). This problem was tackled tentatively in references [TT1 IT2*] in 
what is termed the hydrostatic HNC approach and this will be commented on in section 

ED 

On the other hand, one is inclined to link this failure to the difficulty of 
incorporating the wall as a second component. As the formal framework embodied in the 
graphical expansions holds equally well for a general inhomogeneous situation, one may 
abandon the mixture idea and study the one-component fluid using the common closure 
relations but now for the inhomogeneous correlation functions. These are defined in the 
presence of an external potential which necessarily entails a loss of symmetry for the 
correlation functions and thus results in a noticeable increase in computational power 
required for a numerical solution. Nevertheless, the problem of defining the chemical 
potential in the inhomogeneous situation still persists; for a review of possible solutions 
for a range of commonly used closures see reference Away from phase transitions, 
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inhomogeneous integral equation closures may give very accurate density profiles; for 
the example of a fluid confined in a slit see reference [Tl]. However, systematic studies 
of wetting and drying phenomena within this approach have remained fragmentary, see 
references [T5J HH1 UZ1 UE] • Very recently and quite distinct from the usual closures, an 
implementation of HRT for the inhomogeneous problem has been presented in reference 

Eg. 

Parallel to the development of integral equation approaches, density functional 
models have become established as a widespread tool for the investigation of 
inhomogeneous fluids and for the study of wetting phenomena, in particular. The first 
study of the latter using a 'modern' density functional theory has been reported in 
reference [20]; for a summary of the initial development of the field see reference In 
contrast to integral equations, a very simple mean-field type of functional (with a local 
and therefore actually quite inadequate treatment of the sharply repulsive interatomic 
potentials) predicts rich wetting phenomena at walls, see e.g. references [221 EH] • This 
is related to the underlying bulk equation of state which can be derived from the mean- 
field functional: it is a reasonable zeroth-order approximation to the true equation of 
state. The drawback of the simplest models of the type discussed in references |22J E3] 
consists in their complete failure to describe realistically fluid correlation functions. 
This is due to the inadequate treatment of the short-range repulsive forces. By now 
the treatment of short-ranged correlations within hard sphere density functionals has 
reached a very satisfactory level of precision and internal consistency; for landmark 
developments see references (2U 123 EE] • By contrast, the treatment of the attractive 
tails in the interatomic potentials has remained on the mean-field level which leaves the 
equation of state treated in zeroth order approximation. 

In order to relate density functional results in the presence of a wall, say, to 
quasi-exact results as obtained from simulations, a number of methods have been 
proposed. First, one can observe that away from the critical region, density functional 
results accurately reproduce simulation results if the mean-field coexistence curve is 
appropriately scaled |2Z|- Second, close to exact correlation functions obtained from 
integral equation methods may be used as input in density functionals, see references 
|28| 129} 130] . These amendments generally improve the agreement between density 
profiles obtained from theory and simulation. However, by introducing these ad-hoc 
modifications some key properties of the basic functionals (hard sphere functional plus 
mean field treatment of attractions) are destroyed such as e.g. the compliance with the 
results of exact sum rules. Additionally, if this modified class of density functionals 
is applied to reproduce the bulk fluid correlation functions (this so-called test particle 
limit is obtained by choosing the fluid interatomic potential as the external potential) 
there is, in general, no agreement between the integral equation input and the density 
functional output. 

This requirement of test particle consistency was the starting point for the density 
functional theory developed in [HI] which tries to combine standard density functional 
and integral equation approaches more consistently. By making a Taylor expansion of 
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the free energy functional around bulk densities the standard bulk integral equations 
are derived and the bridge function in the bulk is determined via a density functional 
for a suitably chosen reference system of hard spheres, thereby providing the necessary 
closure relation. The latter resembles the RHNC closure. Results for the bulk structure 
of the one-component plasma and for a few state points of the Lennard- Jones fluid 
[33] show a similar level of agreement (if not better) with simulation data as RHNC 
results. 

In this paper we will develop this idea, the reference functional approach, in 
more detail, generalizing the method to an arbitrary inhomogeneous situation and also 
analyzing bulk fluid and wall-fluid correlations using the mixture description. We will 
see that by introducing the reference functional insertion free energies (i.e., the chemical 
potential in the bulk) are obtained quite naturally. For the wall-fluid problem treated 
in the bulk mixture approach this implies that the insertion free energy of the wall is 
explicitly given; using it we can show that in the case of a hard wall complete drying is 
predicted. 

The paper is structured as follows: In section El we consider a fluid mixture in 
an arbitrary inhomogeneous equilibrium state and derive a general closure for the 
corresponding integral equations using the concept of the reference functional. This 
permits us to derive expressions for the insertion free energies corresponding to the 
change in grand potential when inserting the mixture particles into the inhomogeneous 
equilibrium state. In section |3J we apply the approach to a mixture of a Lennard- Jones 
fluid and hard particles. In the limit of infinite radius of the hard particles (upon 
which the particles become hard walls) and infinitely small hard particle density we 
solve the equations for the fluid and wall-fluid correlation functions. For the bulk fluid, 
we find very good agreement with simulation for both the virial and internal energy 
equation of state. The liquid-vapour coexistence curve is obtained by (i) the usual 
Maxwell construction and [ii] the requirement of complete drying at the hard wall. 
We discuss the form of the wall-fluid bridge function required to account for complete 
drying. Furthermore, we consider the consistency of our theory by analyzing two sum 
rules appropriate to the density distributions at a hard wall. We discuss the relation of 
the reference functional method to commonly used density functional mean field models. 
Section 0] contains a summary and perspectives for further research. 

2. Theory 

We consider a liquid mixture which contains n components. The interaction potential 
between atoms of species i and j is given by u l ^(r) with r = |r» — rj\ denoting 
the distance between the two atoms (i,j = l...n). The 2 th element of the vector 
p(r) = {pi(r), . . . , p„(r)} describes the density distribution of species i and likewise 
the i th element of the vector V(r) = {Vi(r), . . . , V n (r)} denotes the external potential 
acting on an atom of species i. In equilibrium, there exists a unique correspondence 
between the external potential V(r) and the density distribution which we denote by p v 
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to highlight this correspondence. As a consequence, there exists a unique free-energy 
functional [33] 

T[ P ] =r d [ P } + ^[p} (i) 

(3F d [p\ = Wrfr p,(r) (log(p J (r)A J 3 ) - l) (2) 

i J 

which is usually split into the exactly known ideal part (containing the species-specific 
de-Broglie wavelength Aj) and the excess contribution. The grand free energy functional 
is defined by 

fl[p] = T[p\ -Jdr (p id ( Pl , ) + p cx (p,,o) - V(r)) • p(r) (3) 

with p ld (pifl) denoting the chemical potential for an ideal mixture and p ex (p ii0 ) its excess 
(over ideal) for the asymptotic densities p = {pi,o, • • • ,Pn,o}- The asymptotic densities 
correspond to the homogeneous density distribution for V = 0. The components of the 
ideal chemical potential are given by f3p ld = log(pj i0 Af) and (3 = l/{k B T) is the inverse 
temperature with Boltzmann's constant k B - The equilibrium density distribution follows 
from 

5p(r 

The free energy functional generates the hierarchy of (inhomogeneous) direct 



= . (4) 

P=Pv 



correlation functions by 







= -c^-^(v h ...,v n ;p v ). (5) 

P=Pv 



The first member of the hierarchy is determined through equation (jlj): 



c^r; p v ) = log - 0pf (p ) + 0V 4 (r) . (6) 

Po,i 

We are now interested in the change of the density distribution p v generated through 
a perturbation corresponding to the external potential V = V + AV. To this end we 
perform a functional Taylor expansion of the unknown excess free energy in the function 
variable Ap(r) = p(r) — Pv( r ) an d define: 

Here, J 7 ™ (the meaning of the label will become clear shortly) contains terms up to 
second order in Ap and all terms beyond second order are subsumed in JF B . By virtue 
of definition (JHJ), 

/3^ HNC = /M cx (p v ) -J2[ dr c(1) '>; Pv)AMr) - 

\Y.j dv dr ' c(2) ^( r ' r > Pv)Api(r)Ap,(r') , (8) 
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where A ex (p v ) is the excess free energy of the equilibrium state corresponding to the 
external potential V. Then the perturbed density distribution p v , corresponding to the 
external potential V + AV is according to equations (jlj) and (0) 



log E^i^ + 0AVi = [ dv'c^\v,v'-p y )A Pk {v' 



k 

5F B 



(9) 

P=Pv' 



It may appear that we have not gained very much from this formal manipulation: 
the density profile p v , is just expressed through the unknown density profile p v , the 
corresponding unknown inhomogeneous direct correlation function c( 2 )' y '(r, r'; p v ) and 
the unknown functional T B . Therefore we go to the test particle limit, i.e. we define the 
perturbation to V by fixing one atom of species j at position r , thus AV^(r) = u^{r— r ). 
For this particular choice of the perturbing potential, the density profiles are connected 
to the inhomogeneous pair correlation functions and = g l i — 1 through 

Pv,i(r)| y/(r)=v/i(r)+Ml3(r _ ro) = p v ,i(r) g ij (r, r ; p v ) , (10) 

A PV'A r )\v:(r)=V 1 (r)+u^(r-r ) = PV,i( T ) ^( r > r 0! Pv) ■ (H) 

The inhomogeneous pair correlation functions are linked to the direct correlation 
functions via the inhomogeneous Ornstein-Zernike relation 

^(r,r ;pv)-c (a ^(r,ro;p v )=E / dr ' Pv, fc (r')c (2) ^(r, r'; p v ) • 

k 

^(r',r ;p v ). (12) 

Inserting equations (jl(Jj) - (jl2|) into equation Q we recover the general closure relation 
for the inhomogeneous correlation functions 

log^'(r, r ; p v ) + Pu ij (r - r„) = h ij {r, r ; p v ) - c (2) ' ij (r, r ; p v ) - 

5F B 







(13) 

Pi=pv,i 9 lJ 



which permits us to identify the density derivative of T B with the inhomogeneous bridge 
function: 



6«(r J r ;p v ) = /9 



(14) 

Pi( r )=Pv,i(r) g ! J(r,r ;p v ) 



So far all is exact, and at this point we recognize that upon setting T B = we recover 
the inhomogeneous HNC equations. Thus the HNC equations are generated by the 
second-order functional for the excess free energy, JF HNC , in the test particle limit. All 
higher orders in the density expansion of jF ex contribute to the bridge function. 

The unknown density profile Pv( r ) is connected to the correlation functions through 
the exact YBG equations: 

mogpv,i(r)=PVVi(T)+PY, J dr ' /^(r,r';p v )p Vj (r') W^r') .(15) 
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Using the inhomogeneous Ornstein-Zernike relation, equation (|12j) . the last equation 
can be transformed into one connecting pv,« and c^'^: 

Vlogp Vli (r) =(3VV t (r) + Y^ [ dr' c^(r, r'; p v )Vp Vji (r / ) . (16) 

Alternatively, this equation follows directly by taking the gradient in equation (jHJ). Thus 
we see that upon specification of a suitable model for JF B , we have a closed system of 
equations for h l \ c^'^ and py,i consisting of equations (|T2|) . (fT3|) . and (fT5j) or equations 
((121), ©I an d ([IS!) - This set of equations is standard in the theory of classical liquids 
and is usuallly derived using diagrammatic expansions [2] or functional methods [3]. 

The excess free energy pertaining to the external potential V, A ex (p v ), and the 
corresponding one-body direct correlation function c^' l (r; p v ) do not appear in this 
system of equations. The latter can be determined using the Potential Distribution 
Theorem [35] which states that — cW ,l //3 is equivalent to the insertion free energy, i.e. 
the chemical potential, of a particle of species i into the inhomogeneous system defined 
by V. By definition, this insertion free energy is the difference in grand potential with 
the test particle fixed at r , say, and the grand potential without the test particle: 

AufXro; p v ) = (3 fi[p]| p=Pv , - (3 ^[p]| p=Pv = -c (1) <*(r ; p v ) (17) 

where V-{r) = Vj(r) + u^(r — r ). Explicit evaluation yields (using the definition of Q, 
equation (j3J), the expressions for the excess free energy, equations (J7|) and (jBJ), and the 
definitions for the inhomogeneous correlation functions in equations (fTTH) and (|TT|) ) 



^ r °) = Yl / dr i° v jW r ) (log^ l (r, r„) + (3u ji {r - r )) - 
j 

\Y.j dY j dr ' PvjtfPvA*') c (2)jfc (r,r')^(r,r )^(r',r ) 



jk 

I dr p v>j (r) h"(r, r ) + (3 ^[ P ) ]pj{r)=pvd gji{r , ro) 



3 

For compactness, we have suppressed the variable "p v " m all correlation functions 
which indicates that these pertain to the inhomogeneous density distribution defined by 
V. We note that in deriving the last equation, A ex and c- l >' 1 have dropped out on the 
right hand side. To simplify this expression, we use the general closure relation (fT3j) . 
the inhomogeneous Ornstein-Zernike equation (|12j) and the definition (J14I) : 



- c«'(r ) = - cgfc(r ) - £ / dr f»A r ) *>) &*(r, r ) 

j J 

-B 



[P}\ O Hi ^ • ( 19 ) 

L ^ J lpj(r)=pv,j 9 J (r,r ) v ' 

This result, which is the main formal result of our paper, constitutes the generalization 
of the well-known HNC result for the chemical potential in inhomogeneous systems || to 

|| See reference for a derivation using the well-known charging method. 
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closures with non-vanishing bridge functions. Here, — c^c( r o) / P describes the position- 
dependent HNC insertion free energy for J-" B = 0: 

-effic(ro) = E/ rfr Pv.W (~^(r,r ) [**(r, r ) - ^(r, r )] - c^r, r„)) . 

j 

(20) 

Note that c (1)ii is determined by the density distribution p v through equation © which 
imposes a consistency constraint on all subsequent approximations to JF B . Considering 
the limit |r | — > oo, we may assume that the external potential vanishes, V(r ) — > 0, 
thus p v — > p and the inhomogeneous correlation functions become the bulk correlation 
functions, g^ l (r, r ) — > <7 Jl ( r — r o)- Then equations (|BJ) and (fTTIJ) imply an equation for 
the excess chemical potential of species i in the bulk. 

/wao =^r HNC (Po)-E^' / ^^( r )^( r )+/5^ B [p]| Pj(r)=P0j9 , (r) , (2i) 
^r HNC (p ) = / rfr (^'w - c(2)Ji ( r )] - c(2)J ») • ( 22 ) 

The same equation is derived if V = is assumed from the outset, i.e. if the functional 
expansion of the excess free energy is performed around the bulk state described by p . 

In general, the excess free energy functional beyond second order, the bridge 
functional JF B , is not known. However, we may approximate J-" B by a density functional 
for a reference system in the following manner: 

F B [p] « F B ' rcf [p] = F rcf [p] " ^ HNC ' ref [p] , (23) 
where the second order HNC contribution is given by equation (jSJ): 

^HNCref [p] = ^"^(py) " £ / * ^''^ Pv)^P,(r) - 

i E / rfr rfr ' c 2'>> r '; p v )A Pi (r)Ap i (r') , (24) 

For an explicitly given functional jF ref , the direct correlation functions c£f and c^' 2 "' 
are obtained through equation (jSJ). The reference functional approximation implies that 
all direct correlations beyond second order are well approximated by the corresponding 
direct correlations of the reference system. Equivalently, it states that the contributions 
to the excess free energy beyond those of the reference system are of second order in the 
density difference Ap, 

AA = ^ cx [p] - ^ cx - rcf [p] = ^ HNC [p] - .F HNC ' rcf [p] (25) 

= const. - i E / dr (c«»* - cg^r; Pv)Ap,(r) - 



h E / dr dr ' ^ - c S' iJ )(r, r'; p v ) A Pi (r)Ap,(r') . 



(26) 
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Thus the reference functional approximation defines a kind of a generalized mean-field 
approximation in treating a general fluid beyond a given reference system where the 
coefficient in the quadratic term, determined self-consistently. 

In practice, there are only three reference systems for which we possess reasonably 
accurate density functionals. They are all related to hard-sphere systems where we 
can use geometric arguments for their construction. These encompass (i) functionals of 
Rosenfeld type for hard-sphere mixtures (see e.g. references fI5\ l2o]). (ii) a functional 
for the Asakura-Oosawa colloid-polymer mixture [HE], and (Hi) a recently proposed 
functional for non-additive hard-sphere mixtures [37] . The functionals of type (i) have 
been tested in a variety of inhomogeneous situations and can be regarded as very robust. 
The functional (ii) has been applied to the structure of free interfaces and near a hard 
wall, predicting a wealth of interesting phenomena [HE]. A drawback is that the polymers 
are treated only linearly in their density, consequently correlations from the functional 
are deficient if they are sensitive to higher orders in the polymer density. The functional 
(Hi) has not been tested yet for inhomogeneous situations. It contains the functional 
for the Asakura-Oosawa model as a limiting case, so similar merits and problems may 
be expected there. 

Summarizing the literature which is concerned with the reference functional idea, 
we note that the formalism presented above has been formulated previously only for 
functionals of type (i) and for functional expansions around bulk densities p , i.e. 
V = 0. In fact, the reference functional approximation as expressed by equation 
was introduced in reference (HI] and used to describe the one-component plasma with 
high accuracy and its similarity to the reference HNC equations was also pointed out. 
Later on, the numerical efficacy of this idea was illustrated for a few state points of the 
one-component Lennard- Jones (LJ) fluid [33] and of binary mixtures with their size 
asymmetry being not too large [HE]- Comparison with simulation results for both the 
correlation functions and thermodynamic properties shows very good agreement. 

Although the hard sphere functionals of type (i) are very robust, they are not 
exact. For example, the bulk direct correlation function obtained from the functionals 
by direct functional differentiation, equation (jSJ), represent good approximations for its 
behaviour inside the hard core but predict it to be zero outside. In order to improve 
upon the accuracy of the correlation functions, the reference functional method may 
be applied to the hard sphere reference system itself. As a result, for one-component 
hard spheres the bridge function for the bulk integral equations generated by using the 
White Bear functional of reference appears to be the most accurate one known 
up to now. Applied to binary hard-sphere mixtures and size asymmetries > 10, the 
bulk integral equation closure shows deficiencies in calculating the distribution function 
between the big spheres [3JJ]- This may be related to the inadequacy of the functional 
Taylor expansion around bulk densities for large size asymmetries or to the inadequacy 
of the functionals themselves (as manifest in a low-density analysis of the third-order 
direct correlation function |40j). 

The advantage of the present formalism is that for a general, inhomogeneous density 
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distribution p v it is possible to determine directly insertion free energies using equation 
(I19J1 . Of much interest is e.g. obtaining the insertion free energy for colloids at fluid 
interfaces which determines their adsorption behaviour at the interface. However, a 
numerical solution of equations ()12j ) -([15 p requires large numerical efforts. Functional 
minimization in two or three dimensions is complicated by the nature of the three types 
of functionals mentioned above; these contain 5 and 6 1 like distributions. 

In this paper, we confine ourselves to an illustration of the approach for a simple 
one-component fluid and for V = 0. On the one hand, the solution to the homogeneous 
Ornstein-Zernike equation with the integral equation closure defined by equation (fT3j) 
gives the pressure p(p) and the chemical potential p,(p) for a fixed temperature via the 
virial equation of state and equation (|21j). respectively. This permits the determination 
of the liquid-vapour coexistence curve via a standard Maxwell construction as is usually 
done in integral equation theories. On the other hand, by considering a binary mixture of 
the fluid with hard spherical cavities in the limit of vanishing cavity density and infinite 
cavity radius we can derive a wall-particle integral equation. The onset of complete 
drying at the planar hard wall (the cavity with infinite radius) permits an independent 
determination of the coexistence curve. We examine the degree of consistency between 
these two different routes. 



3. Application: Equation of state for a simple fluid and drying at a hard 
wall 

We consider a binary mixture of cut-off and shifted Lennard- Jones particles (species 1) 
with hard cavities of radius R c (species 2). Let 

u LJ (r)=4e[(r/a)- 6 -(r/a)- 12 } (27) 

then the interaction potentials in the mixture are given by 

„n (r) _j«u(r)-^(rj (^j.pg) 

exp(-/V 2 (r)) = exp(-/3w 21 (r)) = 6(R C - r) , (29) 
exp(-(3u 22 (r)) = 6(2R C - r) . (30) 

We apply the formalism from the last chapter with p v = {p^o, P2,o} an d for the reference 
system we use the original Rosenfeld functional which contains two unknown hard- 
sphere diameters d\ and d 2 . In the limit p2,o — * (pi ; o = Po) the equations for h 11 = h 
and c^' n = (fluid correlations) and h 12 (cavity-fluid correlation) decouple. The 
integral equation system for the LJ fluid (cf. equations fll2(l and Q13j) ) is given by: 



B,ref 



(31) 

p(r)=p s(r) 



logp(r) + (3u n (r) = h{r) - c^(r) - (3 § ^ 

h(r) -c (2) (r) =p h*c {2 \r) (32) 

\h * c (2) (r) = J rfr'c (2) (r - r') h(r')j , (33) 



Integral equations for simple fluids in a general reference functional approach 



11 



with g(r) = g{r) etc. Following the reference functional assumption, T B has been 
replaced by the corresponding functional for the reference system, cf. equation 
and the fluid bridge function is given by 



X -77B,ref 

6(r) = r 



(34) 

p(r)=p g(r) 



5p(r 

Using the Ornstein-Zernike equation in the limit p2,o — » in equation ()13|) yields a 
single equation for the cavity-particle correlation which, however, needs as input the 
solution for c- 2 \ the pair direct fluid correlation function: 



log# 12 (r) + (3u 12 {r) = p h 12 * c (2) (r) - (3- 



(35) 

p(r)=P0 S 12 (r) 



5p(v) 

In taking the planar hard wall limit R c — > oo, we shift the origin of the coordinate system 
such that lim^^oo exp(— f3u l2 (v)) = 9{z). In this limit, the cavity-fluid correlation 
functions will be referred to as the wall-fluid correlation functions. 

Note that in the dilute limit p 2 ,o — * both the system of equations (|31|). 
and equation ()35|) can be derived directly from minimization of the grand potential 
employing the excess free energy functional of the one-component L J fluid, ^ rHNC [p(r)] + 
jF B ' ref [p(r); di], and using u 11 and u 12 , respectively, as external potentials. Consequently, 
the reference functional ^ B ' ref is derived from the one-component Rosenfeld functional 
and contains only one unknown parameter, the hard-sphere diameter d±. 

Both equations (|31|) and (|3*3j) are similar to the reference HNC closure jSj. There, the 
bridge function is approximated directly by the bridge function of the reference system 
whereas here the bridge function is obtained from the reference system free energy 
functional. For state points near the triple point, both prescriptions agree with each 
other closely. However, in reference HNC the fluid bridge function remains short-ranged 
near the bulk critical point as does the wall-fluid bridge function near coexistence on 
the liquid side. Thus reference HNC omits effects of bulk criticality and of drying at the 
hard wall. In the present scheme, the fluid bridge function becomes long-ranged through 
its functional dependence on the fluid correlation function g (which is long-ranged near 
the critical point) and so does the wall-fluid bridge function near coexistence through 
its dependence on the wall-particle correlation function g 12 . As will be seen, this does 
not lead to the correct bulk critical behaviour for the fluid but drying at the hard wall 
is described very well. 

For the determination of the unknown hard-sphere diameter di we propose the 
following criterion, 

^ (^ B ' rcf [Po g(r);d 1 } - ^ [p g Icl (r); d,}) = , (36) 

which corresponds to extremizing the free energy difference between the fluid and the 
reference system with respect to d\. The correlation function of the reference system 
g rcf is taken to be the Percus-Yevick correlation function consistent with the Rosenfeld 
functional. This criterion is local in the bulk state space (p, T) and turns out to be 
practicable and reliable, furthermore it is seen in the numerical results that in the 
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moderate to high density region it is equivalent to the usual reference HNC criterion [S]. 
However, a deeper understanding of equation (jBTIj) with regard to a unique specification 
of the reference system is absent at this point. 

For the numerical calculations, we use the cut-off r c = 4 a for which an extensive 
body of simulation data for the equation of state exists |41j . 

3.1. Equation of state 

The pressure p and the excess (over ideal) internal energy per particle u are determined 
from the fluid pair correlation function through the well-known equations: 

— = l--7rpo/3/ drr 6 g(r) — , (37) 

Po 3 J dr 



Figure 1. (colour online) Left panel: Virial equation of state at various temperatures, 
p* = p*(pq). The pressure is calculated according to equation (|57J) . Right panel: 
Internal energy equation of state at various temperatures, (3u — (3u(pq). The excess 
internal energy is calculated according to equation l|38JI ■ The symbols denote simulation 
data taken from reference 01] . 
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POO 

u = 2n p / dr r 2 g(r) u u (r) . (38) 



o 



The usual consistency check of integral equation closures proceeds via the compressibility 
relation 

pdp =1 _ An r dr r 2 c ( 2) ^ _ (39) 







dpo 

In principle, the excess Helmholtz free energy per particle a ex can be determined through 
integration of equation (|3~%j) along isochores, 



pa ex = I u{(3')d(3\ (40) 



o 



and the excess chemical potential follows by differentiating on an isotherm 

M ex = a ex + poj— ■ (41) 
opo 

However, in the reference functional approach the excess chemical potential is given 
directly by (see equation (f2"TJ0 

POO 

(3p ex =/?// x ' HNC -47T Po / drr 2 b(r)g(r) + PT B ' rei [p g(r);d 1 ] ,(42) 

Jo 

^ex.HNc = 4?r pQ J°° dr r 2 (±h(r) [h(r) - c (2) (r)] - c {2) (r) J . (43) 

Consistency with the virial equation (J3*7j) may be checked using the isothermal Gibbs- 
Duhem relation 

^ = Po , 0* = /i id + /^ ex ) • (44) 



From the simulations of reference j4T] the critical temperature and density of the 
cut-off and shifted LJ fluid can be estimated as T* w 1.25 and p* ~ 0.31^[. In Figure ^ 
we present our results for p*(po) and 0u(pq) for temperatures T* = 0.8 . . . 1.4. We see 
that for all temperatures and densities (both close to and away from the critical region) 
the pressure and internal energy calculated from the reference functional match the 
simulation data very well. This indicates that the pair correlation function g{r) is very 
accurate within the range of w n (r) (see equations (|3T|) and (J3~H|0 . In Figure |2l we check 
the consistency of the virial and the compressibility equation of state (upper panel) and 
the Gibbs-Duhem relation (lower panel) with the chemical potential calculated via the 
insertion method, equation (J42|) . Clearly, for the critical region, roughly delineated 
by T* = 1.2 . . . 1.3 and pg = 0.15... 0.5, we observe inconsistencies between the 
virial and compressibility equation of state, and, somewhat weaker, deviations from 
(dp /dp) / po = 1. However, outside the critical region and for stable bulk densities po, 
the ratio of dp/dp calculated via equation (J3l?j) and via the derivative of equation (|3*7j) 
indicates violations of at most 20%, and the Gibbs-Duhem relation is satisfied rather 
well. This will be important for the subsequent investigation of drying at the hard wall. 

If Reduced units, denoted by an asterisk, are defined by setting e = a = 1. 
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Figure 2. (colour online) Testing the consistency of the equation of state at various 
temperatures. Upper panel: The ratio of the density derivatives of the pressure, 
calculated (i) using the compressibility equation (|39J) and (ii) via the derivative of 
a polynomial fit to the virial pressure, equation (|37J) . Lower panel: The ratio of the 
chemical potential derivative of the virial pressure and the bulk density. In order 
to perform the derivative, both the virial pressure, equation H37(l . and the chemical 
potential, equation Ij42(l . were fitted to polynomials in the bulk density. In both 
panels, the density coordinate of the filled circles for the three lowest temperatures 
T* indicates the corresponding coexistence density of the liquid phase obtained from 
simulations |41|. 
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3.2. Two routes to liquid-vapour coexistence and drying at the hard wall 

Integral equations usually give results for pair correlation functions in the metastable 
as well as in the stable region, and therefore it is customary to obtain the coexistence 
curve by requiring p(pi) = p(p g ) and p(pi) = p(p g ); this is equivalent to employing a 
common tangent or Maxwell construction. The density of the coexisting liquid is given 
by pi, while p g denotes the coexisting gas density. The pressure is given by equation 
(jHZj) and the chemical potential may be calculated by either equation or in our case 
equation (}4*2j) . 

However, one may avoid the problem of dealing with correlation functions in the 
metastable domain altogether if one defines p\ as the bulk density for which complete 
drying at a hard wall occurs at a given T < T c . Correspondingly, p g is then the density 
in the bulk of the infinitely extended drying film. 

In order to understand the physics of drying let us consider a liquid of density 
Po near a hard wall at fixed temperature which is close to coexistence, expressed by 
5p = p(po) — p(pi) > and 5p small. Since the wall exerts no attractive forces on the 
fluid molecules, we would expect the fluid density near the wall to be gaslike as this 
is energetically more favourable. Corroborating this argument, there exists a sum rule 
relating the fluid density at the hard wall, p c , to the bulk pressure of the liquid, 



(see also equation below). Since the second virial coefficient of a simple liquid, 
is negative, this sum rule tells us that p c < p g for 5p — > 0. Therefore the density profile 
must pass from a value smaller than the coexisting gas density at the wall to the liquid 
density in the bulk. Since we can safely assume (for T < T c ) that the transition from 
p g to pi happens within a few molecular diameters we are led to the hypothesis that 
upon 5p —>■ a gas layer forms between the hard wall and the bulk liquid whose width 
I goes to infinity as 5p — > for all T < T c - complete drying occurs. More formally, the 
emergence of the gas film may be treated using a coarse-grained interface Hamiltonian 
which in mean-field approximation predicts a slow divergence of the gas film width 
upon approaching coexistence: / oc — In 5p. For a recent general analysis of liquids with 
short-ranged attractions near spherical hard walls, see reference |4*H] . 

The coexisting densities in the hard wall route are obtained by inspecting the grand 
potential 



with u 12 (r) — > u 12 (z), the planar hard wall potential. The excess free energy functional 

^HNC + jrBM 

is given by the Taylor expansion around the liquid bulk density po = p\. 
When examining the stability conditions of an infinitely extended drying film between 



Pc = Pp(Kp)) 

~ /3(p(p(p\))+ Pi5p) , 
~ Pg + P(Pif£/2 + PiSfi) , 



(45) 



Q[p] = ^ id [p(r)] + ^ HNC [p(r)] + ^ [p(r)] 




(46) 
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the wall and the infinitely extended bulk liquid, we can neglect the hard wall potential 
u 12 and the grand potential functional in equation (|4*U|) simply becomes the bulk grand 
potential in a Taylor expansion around p\. As a first condition, the grand potential of 
the drying film with density p g must be stationary 

Mr) pW=p , 

Let V denote the system volume. We introduce the bulk densities of the grand potential, 
cu(p) = Q(p)/V and of the free energy, f(p) = T{p)jV . The second condition is that 
the drying film can only coexist with the liquid if the grand potential densities for the 
two bulk densities p g and p\ are equal: 

u(p g ) - u(jn) = . (48) 

We insert the grand potential of equation (|4~E|) into the above two conditions, evaluate 
them at the required coexisting bulk densities and using the abbreviations A/i = 
p(p g ) — p(pi), AdJ = ui(pg) — uj(p{) we obtain: 

= fJAp = /3p ex ' ref (p g ; d n ) - (3p ex ^( Pl ; d ltl ) + 

[c( 2 \0; Pl )-c^(0;d hl )} (p g — pi) + log — , (49) 

Pg 

= Au = r x < ref ( Pl ; d u ) - / ex ' ref (p g ; d ltl ) + ifo - Pg ) - 

[p 1 /i ex ' rcf (p 1 ; C / lil )-p g /i ex ' ref (p g ;rfi,i)] + 

^[^pd-^iO'M] {pI-pD (50) 

Here, d\ t \ = rfi(pi) is the reference hard-sphere diameter which is determined by 
condition (J36|) and cS 2 '(k) is the Fourier transform of c^ 2 \r). The system of equations 
(|4~9*|) and (|3Ujl can be solved for p\ and p g along with the simultaneous solution to the 
integral equations (p?T]) and for the liquid (which is needed to obtain c^ 2 \0; pi)). 

In Figure E3 we compare the coexistence curve from the standard Maxwell 
construction (using equations ()37|) and ([42))) and that from hard-wall drying with a 
curve fitted to simulation data from reference jH] which was obtained from an equation 
of state. In agreement with our findings concerning the accuracy of the Gibbs-Duhem 
equation (cf. Figure (J2J)), the Maxwell construction gives an excellent description of 
the coexistence curve away from the critical region. Closer to the critical point, we 
are not able to find a solution using the Maxwell construction. This can be traced 
back to the violation of the Gibbs-Duhem equation. Hard-wall drying yields coexisting 
liquid densities with similar accuracy but there are clear deviations from simulations 
for the coexisting gas densities are clearly visible. The reason for this lies in the error 
in extrapolating the bulk grand potential density, equation (jSUJ), to low densities, using 
only bulk properties of the reference system and the pair direct correlation function of 
the coexisting liquid state point. Nevertheless, the fair agreement between the results 
from the Maxwell construction and from hard-wall drying is most encouraging, and, to 
our knowledge, has not been achieved before within integral equation approaches. 
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Figure 3. (colour online) Liquid-gas coexistence for the cut-off and shifted LJ 
fluid. The curve is derived from an empirical fit to the free energy of the full LJ 
fluid, using mean-field corrections to account for the cutoff and shift in the potential 
|41j . The squares are the coexistence points derived by the Maxwell construction 
using equations (|37|l and (|42|l . The diamonds are the coexistence points obtained by 
considering complete drying at the hard wall expressed through equations H49JI and 
(|SIJ). 




Upon approaching the coexisting liquid density from above on an isotherm, the 
wall-liquid adsorption density T w i diverges and the wall-liquid surface tension 7 w i 
approaches the sum of the liquid-gas surface tension 7^ and the wall-gas surface tension 
7 wg . A mean field analysis of an effective interface Hamiltonian for this situation gives 

7wi(/i) = 7w g (/i) + 7^ + f Ap 5(i 

r w i(/i) = T wg (p) - £ Ap In (-^-j . (52) 

Here, /1 = /i(po) is the chemical potential in the bulk, far away from the wall (S/j, = 
ji — yu(pi)), Ap = pi — p g is the difference of the coexisting densities and £ is the decay 
length of the correlation function in the gas phase. The effective parameter a is of order 
0(/cbT/£ 3 ). Fluctuations are expected to change the above behaviour only slightly; a 
treatment within the interface Hamiltonian approach |Hj yields a renormalization of 
the bulk correlation length: £—>(! + uj/2) £ with uj = k B T/(Air 7^ £ 2 ). 



In 



Ap 5p 



+ 1 



(51) 
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Figure 4. (colour online) Surface tension (top panel) and adsorption density (bottom 
panel) of the cut-off and shifted LJ fluid at a hard wall at temperature T* = 1.0. 
Symbols are integral equation results and curves correspond to a fit near coexistence 
to the functional form given in equations Q51|l and H52J) . setting 7wg(/-*) to a constant 
and neglecting r wg (//). From the fit we extract the bulk correlation length £ 0.45 a 
and a w 0.89 fe^T/cr 3 . The range of the chemical potential covered by the integral 
equation results corresponds to the bulk densities Pq — 0.686 . . . 0.9. 



1 
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We have analyzed the surface tension and the adsorption density for the isotherm 
T* = 1.0. The adsorption density is given by 

POO 

r w i = p h 12 (z)dz , (53) 
J o 

while the surface tension 7 w j can be derived from the general formula (|17p. Its explicit 
form is given in | Appendix A In FigureEJwe show the surface tension and the adsorption 



density obtained from the reference functional theory together with a fit to the mean- 
field prediction given by equations and We observe that the latter describes 
the integral equation results not only close to the coexisting density p* ~ 0.686 but also 
up to densities p$ = 0.8 (5p* ~ 1.5) where we can observe a crossover of the adsorption 
density from negative to positive values which is equivalent to a maximum in the surface 
tension. This result is consistent with the idea that the wall-fluid correlations are derived 
from a generalized mean field functional; see equation ()26j) and the discussion thereafter. 
Note that results very similar to those in Figure 0] were obtained in one of the pioneering 
studies of hard-wall drying [42. , in simulations for a square-well fluid. 

Next, we compare the predictions for the wall-fluid bridge function b 12 (z) using 
the reference functional with those of the reference HNC equations. This comparison is 
shown in Figure El for a state point very close to coexistence (pj$ = 0.6859) and a state 
point deep in the stable liquid domain (p$ = 0.9). In reference HNC, the bridge function 
is derived purely from the reference system (hard spheres with diameter d\ at the hard 
wall) and this bridge function is approximated very well by 



^rB,ref 

^RHNc(^) = P' 



(54) 

p(r)=p 3 12 > rcf (z) 



Sp(r 

Here, g 12 > vei is the wall-fluid correlation function for hard spheres at the hard wall 
obtained by minimizing the Rosenfeld hard sphere functional in the presence of the hard 
wall potential. For the higher bulk density (lower panel in Figure EJ) b 12 and &rhnc sh° w 
very similar, short-ranged behaviour whereas for the density close to coexistence (upper 
panel in Figure EJ) b 12 acquires a long-ranged component of approximately the same 
extent as the drying film (compare with the plot of g 12 in Figure EJ). Such a long-range 
behaviour is completely absent in &r HNC which stays short-ranged and consequently the 
drying film formation is not captured by the reference HNC equations. 



3.3. Comparison with the hydrostatic HNC approach 

At this point it is worthwhile to compare the present approach with the hydrostatic HNC 
(HHNC) equations of references ^] which is to the author's knowledge the only 
systematic attempt to account for wetting (drying) phenomena within integral equation 
theories. Following reference jTT] , we imagine that an atom of species 2 (originally our 
hard cavities) exerts a very weak and slowly varying potential on an atom of species 1. 
Let the slowly varying correlation functions between species 2 and species 1 be denoted 
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Figure 5. (colour online) Wall-liquid bridge functions from the reference functional 
method (red lines) and reference HNC (blue lines) at temperature T* = 1.0. Upper 
panel: Bulk density Pq — 0.6859 w p\. Lower panel: Bulk density p^ — 0.9. The 
wall-fluid correlation function g 12 from the reference functional method is given by 
the dashed line. For the density close to coexistence, g 12 displays a thick drying film 
whose extent is similar to that of b 12 , i.e. about 5 a. By contrast, &rhnc remains 
short-ranged. 
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by h s and g s = h s + 1. Then the HHNC bridge function is 

&HH N c(r) = W(Po 9s(r)) - /3// x (p ) - P<A(r) ^p-(po) • (55) 

Here, p ex (p) is the excess chemical potential of the fluid (species 1) which is supposed 
to be known exactly. This expression for the bridge function is exact for slowly varying 
potentials (hydrodynamic limit). Furthermore we note the thermodynamic identity 

= _ 4?r r drr 2 c (2)( r s = _ g(2)(0) (56) 

dp Jo 

which follows from equations and (jH)) . In order to apply this approximation for 
the bridge function to potentials which vary sharply over molecular distances (i.e., the 
hard wall potential or the fluid potential w 11 ) the authors of references [TT) IT2*] introduce 
a normalized weight function w(r) which smoothes the correlation functions: 

h s (r) = w*h 12 {r), g B {r) = w*g 12 (r) . (57) 
As previously, the asterisk denotes a convolution. Thus the bridge function is 

&HHNc(r) = /3p ex (po ft(r)) - /9p cx (p ) + Po £ (2) (0) w * h 12 (r) . (58) 
The wall-fluid integral equation based on this bridge function, 

log g 12 + (3u 12 = p h 12 * - ^ 2 HNC (59) 
supports a drying solution at the "true" coexisting densities. To see this, let p = Pu 



then within the thick drying film g 2 = p g /p\ and (3fj, ex (p g s (r)) — (3p ex (p ) = log(pi/p 
consistent with equation (jB^j) . The drawback of this reasoning is that one needs explicitly 
the function p cx (p) (which one would like to actually obtain from the fluid integral 
equation) and the introduction of the weight function is purely phenomenological. 
Furthermore, if one applies this form of the bridge function to the integral equation of 
the fluid itself, the results are clearly inferior to the present method or to reference HNC 
Thus there appears to be no route in HHNC to obtain the excess chemical potential 
/i ex (p) self-consistently. 

On the other hand, the formally exact bridge function is given by 
5F B 



',12 







5p(r) 



P(r)=Po g 12 (z) 

= Pp cx [p g 12 (z)] - Pp ex (p ) + Po h 12 * tF>{z) , (60) 

which follows from equations (J7J) and (JHJ). Its approximation in the reference functional 
formalism is given by 

b 12 « /3// x ' ref [p g 12 (z); d x ] - /3/i ex ' rcf (p ; d x ) + p h 12 * c^ rcf {z; d x ) (61) 

Here the excess chemical potentials p cx [-] and p ex ' [•] are given by 

r •rcx[,ref] 

p cx ^[p(r)] = ^- r (62) 

In the case of the Rosenfeld functional, p ex ' ref [-] involves double convolutions with a 
certain set of geometric weight functions Comparing equations ()58j) and (ffiTj) 
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we note that instead of employing the exact fluid chemical potential /i ex and direct 
correlation function S- 2 \ the respective reference system quantities appear and that the 
phenomenological weight function w is replaced by more complicated expressions related 
to the nature of the reference functional. The comparison illustrates that compared to 
HHNC, the bridge function within the reference functional method takes account of the 
short-range correlation much better because of the use of the very accurate Rosenfeld 
functional. On the other hand, through the introduction of the reference system, the 
hydrostatic limit is not fulfilled by the reference functional method; this is reflected by 
the fact that the coexisting densities derived from hard-wall drying (equations (|49J1 and 
are nor perfectly consistent with those from the Maxwell construction. 



3.4- Sum rules and the relation to mean-field DFT 

For the planar hard wall problem, two important sum rules exist: 

^ = 9 U (0) , (63) 
Po 

-^f = r wl (p) . (64) 

The first is the contact sum rule (see equation (|45jl ) which links the pressure p in the 
bulk to the fluid density at contact with the wall. The second is the Gibbs adsorption 
equation and links the surface tension to the adsorption density. As usual, the derivative 
is performed on an isotherm. 

The Gibbs adsorption rule remains valid for local density functionals and functionals 
which are constructed from weighted densities (WDA). The contact rule remains only 
valid for free energy functionals of WDA type [27] . This has rendered such functionals 
a popular tool for the investigation of inhomogeneous matter. However, one should add 
the proviso that in equation (|6Hjl . the pressure p is the one derived from the free energy 
functional for bulk densities, i.e. from the compressibility equation of state. Simple 
fluids are almost exclusively treated in density functional theory using a mean-field 
approximation for the attractive tails (see below) for which the consistency between 
the virial and compressibility equation of state is not very good. Moreover, one might 
object that agreement with sum rule (|H4*j) does not guarantee the accuracy of 7 w i and 
T w i themselves, as these quantities depend strongly on the quality of the free energy 
functional. 

For hard spheres, these problems are solved to a high degree of accuracy by presently 
available functionals [2E1 113] • The simplest mean-field treatment of fluids in which the 
pair potential u u (r) has an attractive tail is given by 

J^\p{r)\ = ^ cx ' rcf [p(r); d x = const.] + \ J dr p(r) w M * p(r) , (65) 

where we consider a one-component system. Here, the weight function w at t is given 
by the attractive tail of the interparticle potential (for which various recipes exist in 
the literature). If we compare the second term in equation (J65|) with the excess free 
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Figure 6. (colour online) The difference between the direct correlation functions of 
the fluid and of the reference system for various stable liquid states. 




r/ c 

energy beyond the reference system in our approach (see equation (J2I3)), then the 
reference functional approach predicts roughly mean-field behaviour for the density 
deviations from the bulk. Thus we may conclude, if the prefactor of the quadratic term 
in equation (J2E1), (c^ — c^' rei )/P, does not vary appreciably with the bulk density p 
then a description in terms of the functional given in equation (Jo^ will capture the 
essential physics. Checking the numerical results (see Figure EJ), we find that for r > a 
indeed a single function describes this prefactor for all stable liquid state points except 
for the ones in the critical region. Inside the harshly repulsive core (r < a) we notice 
some variation which is, however, small when compared to the magnitude of the direct 
pair correlation function of the reference system. Thus the "mean-field" expansion 
around the bulk density as given in equation 1)26)1 is expected to entail similar physics 
as that in the mean-field functional ()U3)) but the former has the effect of giving a more 
accurate equation of state and very reliable fluid correlation functions. The mean field 
approach ()2Tj) is, in turn, computationally less demanding and has the advantage that 
the sum rules hold analytically (with respect to the compressibility equation of state). 
It is important to note that the very concept of expanding around a fixed bulk 
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Figure 7. (colour online) Testing the consistency of the contact density sum rule, 
equation (|63|l . We have plotted the ratio p* /(pq g 12 (Q)) for various temperatures 
T* = 0.8 . . . 1.4. The sum rule value of this ratio, T*, is shown by a horizontal straight 
line. For the three lowest temperatures, results are plotted for the stable liquid phase 
only. 




density implies that the sum rule proofs which work for weighted density functionals do 
not hold for the Taylor expanded functional |46| . Therefore we have tested both sum 
rules numerically for various temperatures. In Figure [7| we have plotted the ratio of 
the virial pressure to the density at contact, p* /[pog l2 (0)} (which should give T*), for 
temperatures ranging from T* = 0.8 . . . 1.4. The deviations are up to 20%. However, 
despite these deviations these results are encouraging since here the virial (= quasi- 
exact) pressure is used. In Figure El we compare the adsorption densities calculated by 
equation (J53j) with the adsorption density obtained from the surface tension via the sum 
rule (jfi4j) . Again we find very reasonable agreement. 

4. Summary and conclusions 

In this article, we have analyzed the pair correlation functions of an inhomogeneous 
fluid mixture by means of a functional Taylor expansion of the free energy around an 
inhomogeneous density profile. We have derived a general integral equation closure for 
inhomogeneous systems which is based on the introduction of a reference functional for 
the excess free energy beyond second order in the Taylor expansion. A benefit of this 
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Figure 8. (colour online) Testing the consistency of the Gibbs adsorption rule, 
equation (|64|l . We show a comparison of the adsorption density T w i — po f dzh 12 
(symbols) with the same quantity obtained from sum rule H64|) as a function of the 
bulk density for various temperatures T* = 0.8 ... 1.4. To perform the /i derivative of 
the surface tension 7 w i, both /i and 7 w i have been fitted to polynomials in po- The fit 
is of limited accuracy only very close to coexistence. The vertical dashed lines for the 
three lowest temperatures mark the coexisting liquid densities. 





approach is that explicit expressions (see equation (JIHJ)) for the insertion free energies 
of particles into the inhomogeneous equilibrium distribution are obtained, which reduce 
to explicit expressions for the excess chemical potential if the expansion is performed 
around a homogeneous (bulk) state. 

As an application, we have considered the equation of state for a bulk cut-off and 
shifted Lennard- Jones fluid and the structure and adsorption of this fluid at a hard 
wall. Treating this problem within the mixture theory and expanding around the bulk 
state, we find excellent agreement between the virial and internal energy routes to the 
equation of state for the fluid and with simulations. Thermodynamic consistency is 
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good for stable state points outside the critical region. By virtue of the underlying 
hard sphere reference functional, complete drying is predicted for the hard wall-liquid 
interface and the corresponding coexisting densities for the liquid branch are in excellent 
agreement with those from the Maxwell construction. However, there are deviations for 
the coexisting densities on the gas branch. Nevertheless, this reasonable consistency 
between the Maxwell construction and hard wall drying is a novel result within integral 
equation approaches. 

Two sum rules related to the hard wall problem, the contact density rule, equation 
(J?)3J), and the Gibbs adsorption equation (JMj) . have been investigated over a wide range 
of densities and temperatures. Although these sum rules are not satisfied exactly, the 
deviations are reasonably small, especially near coexistence. 

The numerical results for both the fluid-fluid and the wall-fluid correlations suggest 
a generalized mean field behavior outside the critical region. In the reference functional 
approach, the difference between the fluid free energy and the free energy of the reference 
system is expanded up to second order in the density difference (from the bulk density), 
see equation ()26|) . The coefficient of the second-order term, the difference between the 
pair direct correlation functions of the fluid and the reference system, is determined self- 
consistently and turns out to be quite insensitive to variations in temperature and bulk 
density for liquid states, see Figure |H1 Therefore we argue that for physical problems 
related to wetting and drying a description with density functionals composed of an 
accurate reference part and a mean field attractive term should be sufficient to capture 
all the essential features and will be much less demanding in computational effort. 

Related to the efficacy of generalized mean-field treatments of simple fluids, we 
note that a quite different mean-field treatment of the density deviations from the bulk 
was already developed in reference [UJ in which the behaviour of a LJ fluid at a hard 
wall was studied as well. In contrast to the present study, the properties of the hard 
sphere reference system were fixed by considering the first YBG equation, for details see 
reference jUJ. Good agreement between simulated and calculated density profiles was 
obtained for selected supercritical state points. 

We anticipate interesting perspectives for further computations in the problem of 
small colloidal particles at interfaces. The free energy profile for particles dragged 
through the interface or their mutual interactions at the interface are only beginning 
to be understood. These observables can be analyzed within the reference functional 
approach for inhomogeneous distributions; the challenging numerical implementation 
remains a task for future work. 
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Appendix 

Appendix A. Surface tension at the hard wall 

The hard wall-liquid surface tension is defined by the expression 

POO 

7 w i = / dz(u(z) +p) , (A.r 



where lu(z) is the grand potential volume density and p = p(po) is the bulk pressure 
of the liquid with density p Q . In terms of the grand potential density functional, this 
equation is equivalent to 

7 wi = (O[p(r)]| p(r)=posl2(z) - n[p(r)]\ p{r)=poe(z) ) I A , (A.2) 

where A denotes the area of the wall. We use the definition 

%(r)] =-F[p(r)]-|rfr(Mpo)-^ 12 (r))p(r) . (A.3) 



with £ = 1 for the first term in equation (jA.2|) and £ = for the second term. According 
to the discussion below equation concerning the hard wall limit, the functional of 
the free energy is the functional for the one-component LJ fluid, 

TW\ = F d [p(v)] + ^ HNC [p(r)] + ^ [p(r); dj , (A.4) 

with di = di(p ) denoting the reference system hard sphere diameter. Straightforward 
evaluation gives 



wl 



dz Po (V 2 hV 2 - h"{l +/3^ ref (po)) - f [(cf - <?W) * h U] h ^ + 

POO 

PF ex > vei [po 9 12 ] ~P dz f e ^ + J wall , (A.5) 
Jo 

fO 



'wan : / dz /?[/ cx (p ) - / cx ' rcf (Po)] + (A.6) 
) V 

M 12 /3(p ex (Po) - P cx ' rcf (Po)) - f [(cf - cf - f ) * /> 12 ] /> 12 ) . 

Here, / ex (po) and / ex ' ref (po) are the bulk free energy volume densities of the fluid and 
the hard-sphere reference system, respectively. We note that the non-local expansion 
of the free energy in terms of Ap = p h 12 (= — p for z < 0) introduces contributions 
inside the wall, contained in J wa ii. Furthermore, 

POO 

cf )[ ' vci] (z) = 2tt / rdrc (2)[ ' ref] (r) , (A.7) 

J\z\ 

and we have introduced the one-dimensional convolution 

POO 

c (2)[,rcf] ^ h 12 ^ = / dz r c (2)[,re{]( z _ ^g) 
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The unknown free energy density of the liquid / ex (po) appearing in the integral 7 wa n is 
fixed by the requirement that beyond the wall (z < 0) the density is zero and therefore 
also / ex (0) = 0. On the other hand, / ex (0) is determined by the functional Taylor 
expansion of the fluid excess free energy around po , evaluated at zero density. Therefore 
we find: 

r(po) = r M ( P o) + Po (^(Po) - ^ cx ' rcf (po)) + | [£ (2) (0) - £ (2) ' rcf (0)] . (A.9) 
This expression is inserted into 7 wa u and we obtain the expression: 

/wan = P i f dz ([c< 2 >(0) - c^ cf (0)] - [(cf - c?)> ref ) * h 12 ] h 12 ) . (A.10) 

* J -oo 

This completes the prescription for calculating 
References 

[1] Morita T and Hiroike K 1960 Prog. Theor. Phys 23 1003 

[2] Stell G 1964 in The Equilibrium Theory of Classical Fluids edited by Frisch H L and Lebowitz J 
L (New York:Benjamin) p 11-171 

[3] Percus J K 1964 in The Equilibrium Theory of Classical Fluids edited by Frisch H L and Lebowitz 
J L (New YorkrBenjamin) p 11-33 

[4] Pastore G, Akinlade O, Matthews F, and Badirkhan Z 1998 Phys. Rev. E 57 460 

[5] Parola A and Reatto L 1995 Adv. Phys. 44 211 

[6] Pini D and Stell G 2002 Physica A 306 270 

[7] Vompe A G and Martynov G A 1996 J. Chem. Phys. 106 6095 

[8] Lado F, Foiles S M, and Ashcroft N W 1983 Phys. Rev. A 28 2374 

[9] Hansen J-P and McDonald I 1990 Theory of Simple Liquids (London:Academic Press) 
[10] Evans R, Tarazona P, and Marconi U M B 1983 Mol. Phys. 50 993 
[11] Zhou Y and Stell G 1990 J. Chem. Phys. 92 5533 
[12] Zhou Y and Stell G 1990 J. Chem. Phys. 92 5544 
[13] Kjellander R and Sarman S 1989 J. Chem. Phys. 90 2768 
[14] Kjellander R and Sarman S 1990 Mol. Phys. 70 215 
[15] Nieminen R and Ashcroft N W 1981 Phys. Rev. A 24 560 
[16] McGough R A and Miller M D 1986 Phys. Rev. A 34 457 
[17] Plischke M and Henderson D 1986 J. Chem. Phys. 84 2846 
[18] Bruno E, Caccamo C, and Tarazona P 1987 Phys. Rev. A 35 1210 

[19] Orlandi A, Parola A, and Reatto L 2004 A microscopic approach to critical phenomena at 

interfaces: an application to complete wetting in the Ising model Preprint cond-mat/0409188 
[20] Ebner C and Saam W F 1977 Phys. Rev. Lett. 38 1486 

[21] Evans R 1992 in Fundamentals of Lnhomogeneous Fluids edited by Henderson D (New York:Dekker) 
p 85-175 

[22] Evans R and Tarazona P 1983 Phys. Rev. A 28 1864 
[23] Kroll D M and Meister T F 1985 Phys. Rev. B 31 392 
[24] Tarazona P 1985 Phys. Rev. A 31 2672 
[25] Rosenfeld Y 1989 Phys. Rev. Lett. 63 980 

[26] Roth R, Evans R, Lang A, and Kahl G 2002 J. Phys.: Cond. Matt. 14 12063 

[27] van Swol F and Henderson J R 1989 Phys. Rev. A 40 2567 

[28] Choudhury N and Ghosh S 2001 J. Chem. Phys. 114 8530 

[29] Gillespie D, Nonner W, and Eisenberg R S 2003 Phys. Rev. E 68 031503 

[30] Tang Y and Wu J 2003 J. Chem. Phys. 119 7388 



Integral equations for simple fluids in a general reference functional approach 



29 



Rosenfeld Y 1993 J. Chem. Phys. 98 8126 
Rosenfeld Y 1996 Phys. Rev. E 54 2827 
Rosenfeld Y 1998 Mol. Phys. 94 929 
Evans R 1979 Adv. Phys. 28 143 
Henderson J R 1983 Mol. Phys. 50 741 

Brader J M, Evans R, and Schmidt M 2003 Mol. Phys. 101 3349 
Schmidt M 2004 Rosenfeld functional for non-additive hard spheres Preprint cond-mat/0407013 
Kahl G, Bildstein B, and Rosenfeld Y 1996 Phys. Rev. E 54 5391 
Oettel M 2004 Phys. Rev. E 69 041404 

Cuesta J A, Martinez- Raton Y, and Tarazona P 2002 J. Phys.: Condens. Matter 14 11965 
Johnson J K, Zollweg A, and Gubbins K E 1993 Mol. Phys. 78 591 
Henderson J R and van Swol F 1985 Mol. Phys. 56 1313 

Evans R, Henderson J R, and Roth R 2004 Non-analytic curvature contributions to solvation free 
energies: influence of drying Preprint cond-mat/0410179 accepted in J. Chem. Phys. 

Dietrich S 1988 in Phase Transitions and Critical Phenomena edited by Domb C L and Lebowitz 
J L (LondomAcademic) vol 12 p 1 

Bryk P, Roth R, Mecke K R, and Dietrich S 2003 Phys. Rev. E 68 031602 

Sweatman M B 2000 Mol. Phys. 98 573 

Weeks J D, Katsov K, and Vollmayr K 1998 Phys. Rev. Lett. 81 4400 



